Testing Lorentz Invariance with neutrino burst 
from supernova neutronization 
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Quantum-gravity (QG) effects might generate Lorentz invariance vioiation by the interaction of 
energetic particles with the foamy structure of the space-time. As a consequence, particles may 
not travel at the universal speed of light. We propose to constrain Lorentz invariance violation for 
energetic neutrinos exploiting the v e neutronization burst from the next galactic supernova (SN). 
This prompt signal is expected to produce a sharp peak in the SN v e light curve with a duration 
of ~ 25 ms. However presence of energy-dependent Lorentz invariance violation would significantly 
spread out the time structure of this signal. We find that the detection the SN v e burst from a 
typical galactic explosion at d = 10 kpc in a Mton-class water Cerenkov detector, would be sensitive 
to a quantum-gravity mass scale Mqg ~ 10 12 GeV (2 x 10 5 GeV) for the linear (quadratic) energy 
dependence of Lorentz invariance violation. These limits are rather independent of the neutrino 
mass hierarchy and whether the neutrino velocity is super or subluminal. 

PACS numbers: 14.60.Pq, 97.60.Bw 



I. INTRODUCTION 

Lorentz invariance violation (LIV) arises in many ap- 
proaches of Quantum Gravity (QG) theories, suggesting 
Lorentz symmetry might be violated at very high ener- 
gies (see [l| for a review). In this context, the space-time 
foam due to QG fluctuations might cause energetic parti- 
cles to propagate at speed v different from the velocity of 
light c, which would be approached only by low-energy 
masslcss particles. This LIV effect can be phenomeno- 
logically parametrized as 
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where n = 1, 2 denote linear or quadratic deviation for su- 
pcrluminal (+) or subluminal (-) particles, respectively, 
occurring at a mass scale Mqg- 

Apart from photons [H, neutrinos from astrophysical 
sources can also be useful to constrain LIV effects. In 
case of ~ O(GeV) energy neutrinos from Gamma Ray 
Bursts or Active Galactic Nuclei Q , since one expects to 
detect at most one or two such neutrinos these sources, 
the main attempt to constrain the LIV scale is through 
the energy dependent time of flight delay of neutrinos 
and the corresponding 7 rays. The limit that would 
be placed is Mqg ft 10 17 GeV for the n = 1 case, 
and Mqg si 10 10 GeV for the n = 2 case, respec- 
tively. Unfortunately till now, no high-energy v from 
these astrophysical sources has been detected. However, 
also core-collapse SNe are guaranteed emitters of MeV 
neutrinos, that have been widely discussed in literature 
to put bounds on LIV. In particular, from SN1987A 
data, exploiting the flight delay of a few hours of the 
7 with respect to v, it has been obtained a bound of 
(v - c)/c < 10~ 12 0, corresponding to M QG > 10 9 GeV 
for n = 1, or Mq G > 10 3 GeV for n = 2, taking 
E ~ 10 MeV as typical SN v energy. 



Moreover, since the expected duration of the SN1987A 
v burst (~ 10 s) was in agreement with the observed one, 
this limits LIV effects which would have broadened the 
SN v time structure due to the energy dependency of the 
Lorentz violation. From the absence of an anomalous 
dispersion, it has been found Mqg ^ 2 x 10 10 GeV for 
the n = 1 case, and Mq G > 4 x 10 4 GeV for the n = 2 
case Q. These bounds arc stronger than the previous 
ones obtained from the SN v time of flight because the 
observed time dispersion is a few tens of seconds, and a 
detailed statistical analysis has been performed. More- 
over, these limits are also much more stringent than the 
ones obtained from the v time of flight measurements in 
the long-baseline MINOS @, OPERA Q and ICARUS 
experiments Q giving Mqg ^ 10 5 GeV for the n = 1 
case, and Mqg 10 2 GeV for the n = 2 case. 

One expects that these bounds could be significantly 
improved with the observation of a high- statistics sig- 
nal from the next galactic SN. At this regard, it has 
been estimated that from the time structure of the v 
signal from a SN at d = 10 kpc observable in the Super- 
Kamiokande detector (with fiducial mass 22.5 kton), it 
would be possible to constrain Mqg ^ 2(4) x 10 11 GeV 
for subluminal (superluminal) propagation in the n = 1 
case, and Mqg > 2(4) x 10 5 GeV for subluminal (su- 
perluminal) propagation in the n = 2 case Ij. An even 
stronger bound (up to two orders of magnitude beyond 
previous estimates) could be obtained exploiting varia- 
tions in time on the scale of a few milliseconds found in 
multi-dimensional SN simulations However, probably 
these features could be detected only for a very close-by 
SN (at d < 2 kpc) 0. 

In the following, we propose to explore the effects of 
LIV on the prompt SN is e neutronization burst. This 
signal is a common feature pd| found in all sophisticated 
supernova simulations. Physically the newly formed SN 
shock disintegrate the dense iron core into free neutrons 
and protons. The electron rich environment behind the 
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FIG. 1: Time evolution of the neutrino number flux for a 15 
Mq SN simulation from the Garching group, for v e (continu- 
ous curve), i>e (dashed curve) and v x (dot-dashed curve). 



SN shock triggers rapid electron capture on these free 
protons producing a huge v e flux. Thus a large number 
of i/ e 's are released when the shock breaks the neutri- 
nosphere and the i/ e 's escape freely. This process of delep- 
tonization 'neutronizes' the environment. The v e and 
muon and tau neutrino-(anti)neutrino (denoted with u x ) 
luminosities rise much more slowly compare to the ^ e 's. 
These facts are very well depicted in Figure [TJ where we 
plot the neutrino number flux of different flavors in the 
neutronization phase for a 15 Mq SN simulation from 
the Garching group The signature of the neutron- 
ization peak is nearly independent from the details of SN 
models, like electron capture rates, nuclear equation of 
state and the progenitor mass. Indeed, it can be con- 
sidered as a "standard neutrino candle" in the simula- 
tions of core-collapse supernovae. In fact probing the 
time-spectra of the neutronization peak one can distin- 
guish the case of inverted neutrino mass hierarchy (IH, 
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< 0) from the normal mass hierar- 



> 0) [11], where the peak is suppressed 
for a value of a 1-3 mixing angle $13 as lar ge a s the one 
currently measured by reactor experiments [12| . [HI ] . The 
LIV inducing an energy dependence in the v velocity will 
broaden the time structure to the v e burst as neutrinos 
with different energy will undergo different time delay. 

In the following we characterize this effect. In Sec. HTl 
we describe the SN neutrino flavor conversions during the 
v e burst phase. In Sec. lIIII we recall the calculation of the 
SN neutrino events rate in a Mton class water Chercnkov 
detector. In Sec. IIVI we present the effects of the LIV on 
the v signal during the neutronization burst. Finally, we 
summarize our results in Sec. [V] 

II. NEUTRINO FLAVOR CONVERSIONS 

The SN neutrino signal observed at Earth is processed 
by flavor conversions during the propagation in the stel- 



lar envelope. Since negligible v e and v x fluxes are emitted 
during the neutronization phase, self-induced flavor con- 
versions are absent (l4l |. Thus, only Mikheyev-Smirnov- 
Wolfenstein flavor transitions occur while the neutrinos 
propagate through the stellar envelope [lj|. Given the 
discovery of large $13 [l2l Il3j the SN v flux at Earth 
is straightforward to calculate. Namely, the v e and v e 
fluxes in normal hierarchy (NH) are given by [21 
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while in inverted hierarchy (IH) one gets 

F Ve = sin 2 tf 12 (^° _F°J+F° 
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Here $12 is the 1-2 mixing angle, with sin 2 $12 ~ 0.31 [T7| 
and is the initial flux for the a'th flavor. The fluxes 
of the other flavors at earth can be found from flavor 
conservation, F Ve + 2F Vx = F° e + 2F° x and F Pc + 2F 9=[ = 

III. NEUTRINO EVENT RATE 



A clear detection of the v e neutronization burst re- 
quires large detectors with a sensitivity to the electron 
neutrino flavors. In particular, a Liquid Argon Time pro- 
jection chamber with a mass of 0(100) kton would be 
able to distinguish this signal via charged-current reac- 
tions on Ar nuclei [l8[. Moreover, also Mton class water 
Chercnkov detectors could get a signature for the v e neu- 
tronization burst, through the clastic scattering of v e on 
electrons (is e + e~ — > is e + e~). In the following, for defini- 
tiveness we focus on a future water Chercnkov detector 
with a fiducial mass of 740 kton [lj|, like the Japanese 
Hyper-Kamiokande experiment. For our calculations, we 
take as benchmark the results of the 15 Mq Garching 
simulations shown in Fig. 1 to characterize the origi- 
nal time-dependent SN neutrino fluxes [F v — F u (E,ti)) 
during the early neutrino emission phase. Here E is the 
neutrino energy and ti is the emission time at the source. 

The arrival time of a neutrino of energy E in the detec- 
tor at earth is = U + d/c+ At, where d is the distance 
of the source. The energy dependent time delay due to 
QG is given by 



At = ± 
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where the (+) sign is for subluminal and (-) for superlu- 
minal case. The constant term d/c is dropped for conve- 
nience. Then the event rates of neutrinos at the detector 
from a SN at distance d, is given by 
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FIG. 2: Neutronization events rate per time bin in a 740 kton water Cerenkov detector for a SN at d — 10 kpc. The upper 
two panels describe subluminal LIV whereas the lower panels are for the superluminal scenario. The right panels are for NH 
and the left ones for IH. The Lorentz Invariant (LI) cases are shown in continuous lines. For LIV we consider the linear energy 
dependence, i.e n = 1. The dot-dashed lines describe LIV for Mqg — 5 x 10 12 GeV and the dashed lines describe LIV for 
M QG = 10 12 GeV. The error bars reflect the number of expected neutrino events in each time bin. 



Here tit is the number of target nuclcons, u{E) the 
elastic-scattering cross section for a neutrino energy E 
on the target [llj. In the following we will assume a 
typical SN at d = 10 kpc. Our calculation of the events 
rate closely follows the one presented in ll| to which we 
address the interested reader for further details. 

IV. RESULTS 

In Figure 2 we show the event rate expected during the 
neutronization burst for different scenarios. In particular 
the left panels refer to the NH cases, while the right ones 
to IH. Moreover, the upper panels show the subluminal 
scenarios and the lower ones the superluminal ones. In 
the absence of Lorentz violation (continuous curves) the 
difference between the time structure of the events for 
the two different mass hierarchies reconfirms the results 
presented in [ll[. In particular, since in the NH case 
the signal is dominated by F® [see Eq. ©] the peak is 
suppressed. Instead, for the IH, the flux is dominated by 



F® [see Eq. ((SJ)] and hence the signal has a peaked time 
spectrum. The standard scenario significantly changes 
once the LIV effects are included. We show examples 
with the linear energy dependence case, i.e n = 1 with 
M QG = 10 12 GeV (dashed curves) and 5 x 10 12 GeV 
(dot-dashed curves). LIV tends to flatten and spread out 
the original time structure in the signal. In particular 
in the subluminal case for Mqg = 5 x 10 12 GeV the 
strength of the peak is reduced and shifted to later times 
(by ~ 5 ms), while for Mqg = 10 12 GeV the peak is 
completely washed out and the signal presents a shape 
that monotonically rises with time. In the superluminal 
cases, for Mqg — 5 x 10 12 GeV the peak is shifted to 
earlier times in both hierarchies (by ~ 5 ms). For NH 
with Mqg — 10 12 GeV the time spectrum becomes flat, 
whereas in IH it monotonically decreases. 

We also considered the case of quadratic energy depen- 
dence of LIV in Fig. 3. Absence of LIV features would 
constrain Mqg 2 x 10 5 GeV in both superluminal and 
subluminal cases. 
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FIG. 3: Neutronization event rate per time bin in a 740 kton water Cerenkov detector for a SN at d = 10 kpc. Here the energy 
dependence of LIV is quadratic, i.e n = 2. The left panel refers to NH case, while the right panel to IH. The cases without LIV 
are represented by the continuous curves. In the cases with LIV we take as representative QG mass scale Mqg = 2 x 10 5 GeV. 
The subluminal cases are in dashed curves, while the superluminal in dot-dashed ones. 



V. SUMMARY 

Wc have studied effects of Lorentz invariance violation 
in the neutrino sector on the SN neutronization burst. 
The analysis shows that the Lorentz invariance violation 
would produce a strong suppression of the expected neu- 
tronization peak in both the superluminal and sublumi- 
nal cases. We find a sensitivity to the quantum-gravity 
mass scale of M QG ~ 10 12 GeV (2 x 10 5 GeV) for the 
linear (quadratic) energy dependence of LIV in both the 
neutrino mass hierarchies. 

This method would give a clean bound on quantum- 
gravity mass scale, better and simpler than the present 
limits from SN neutrinos. Indeed these potential lim- 
its come from the different shapes of the neutronization 



peak, which is a rather model independent feature of SN 
simulations. 
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